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(57) Abstract 



MRI signals detected by a phased array of multicoils (1, 2,..Nc) are combined (20) to provide phase contrast image signals with 
improved signal-to-noise ratio. Detected complex image signals are obtained for different magnetic gradients and are combined to obtain 
phase, magnitude, and magnitude weighted image signals for imaging. 
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PHASE-CONTRAST MRI USING PHASED-ARRAY MULTICOTI. 

5 

A BACKGROUND OF THE INVENTION 

This invention relates generally to phase contrast 
magnetic resonance imaging (MRI) , and more particularly the 
invention relates to phase contrast MRI using a multicoil 

10 phased array for magnetic resonance signal detection. 

As described by Pelc et al. ("Phase Contrast Cine 
Magnetic Resonance Imaging" Magnetic Resonance Quarterly 
Vol. 7 No. 4, pp. 229-254, 1991), phase contrast MRI refers to 
a family of MR imaging methods that exploit the fact that 

15 spins that move through magnetic field gradients obtain a 

different phase than static spins, enabling the production of 
images with controlled sensitivity to flow. Among the 
advantages of phase contrast MRI are its adjustable 
sensitivity to flow and, in some cases, the quantitative 

20 nature of the resulting data. In the most common method, two 
data sets are acquired with a different gradient first moment 
in one direction, and the pixel by pixel phase difference a0 
in the resulting images is used to measure the velocity, V, in 
the direction of the gradient change: a^yaM-lV, where aM x is 

25 the change in gradient first moment and y is the gyromagnetic 
ratio. Four measurements can be used to measure all 
components of velocity. Sometimes, for aesthetic reasons 
related to the very high noise in a0 in regions of low signal, 
"magnitude weighted" or "magnitude masked" velocity images 

30 with pixel intensities proportional to Ma0 are produced, where 
M is the magnitude image. 

Phased array multicoils can be used in MRI to 
improve the signal-to-noise ratio (SNR) of images. See, for 
example, Roemer et al. "The NMR Phased Array", Magnetic 

35 Resonance in Medicine, Vol. 16, pp. 192-225, 1990; Hayes et 
al. "Noise Correlations in Data Simultaneously Acquired from 
Multiple Surface Coil Arrays" Magnetic Resonance in Medicine, 
Vol. 16, pp. 181-191, 1990; and Hayes et al. "Volume Imaging 
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with MR Phased Arrays" Magnetic Resonance in Medicine, 
Vol. 18, pp. 309-319, 1991. Basically, the multicoils allow 
large field of view imaging with the SNR of a small surface 
coil. Separate images are acquired and reconstructed from 
5 each of the elements in the multicoil array. These separate 

images are then combined into a single image, with each coil ^ 
dominating the spatial regions where its SNR is the highest. 
Although various algorithms for this combination process have 
been proposed, the most commonly used methods use the 

10 magnitude of the pixel values in the individual images.. For 
example, one very common method forms images as the square 
root of the sum of squares of the individual images. 

The combination of phase contrast techniques and 
phased array coils would be very desirable in order to improve 

15 the SNR of the velocity-sensitive images. However, a simple 
and efficient method for combining the information received 
from the various coils that retains a proportionality with 
velocity is needed. 

20 SUMMARY OF THE INVENTION 

In accordance with the invention, the complex image 
data (Z ± ) from each coil (i) in a multicoil phased array 
including the magnitude and the phase (<p ± ) are employed 

along with specific coil gain factors (W \) to generate phase, 
25 magnitude, and magnitude-weighted images. 

More particularly, let Z lL and Z 2 ± be the complex 
images reconstructed from coil i with sequences 1 and 2 , and 
magnitude image, M if equal 

M ± « (|Z XJ | + \Z 2i \)/2 
30 and phase image, <f> if equal 

p 

<|) i = phase (Z xi Z 2i *) 

where * indicates the complex conjugate and 

phase (Z) = arctan2 (imag (Z ) , real (Z) ) . 
Further, let VI ^ be the gain factor for coil i. 



WO 95/04947 



PCT/US94/08954 



Then, in accordance with the invention, the phase, 
magnitude, and magnitude-weighted image are given as: 

a0 s arctan2 ( imag(C), real (C) ) 
mag weighted image = M a0 

where 



c = £[( z i; z 2i yw>] 



The invention presents a simpler and more optimal 
method in using phase contrast MRI with a phased array 
multicoil than heretofore known. 
10 The invention and objects and features thereof will 

be more readily apparent from the following detailed 
description and appended claims when taken with the drawing. 

BRIEF DESCRIPTION OF THE DRAWINGS 
15 Figs. 1A-1D illustrate the arrangement of MRI 

apparatus and magnetic fields generated therein. 

Fig. 2 is a function block diagram of MR imaging 

apparatus . 

Fig. 3 is a schematic view of an MRI phased array as 
20 used in NMR apparatus. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
Referring now to the drawings, Fig. 1A is a 
perspective view, partially in section, illustrating coil 

25 apparatus in MR imaging system and Figs. IB-ID illustrate 
field gradients which can be produced in the apparatus of 
Fig. 1A. Briefly, the uniform static field B 0 is generated by 
the magnet comprising the coil pair 10. A gradient field 
G(x)is generated by a complex gradient coil set which can be 

30 wound on the cylinder 12. An RF field B x is generated by a 
saddle coil 14. A patient undergoing imaging would be 
positioned along the Z axis within the saddle coil. In 
Fig. IB an X gradient field is shown which is parallel to the 
static field B 0 and varies linearly with distance along the X 
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axis but does not vary with distance along the Y or Z axes. 
Figs. 1C and ID are similar representations of the Y gradient 
and 2 gradients fields, respectively. 

Fig. 2 is a functional block diagram of conventional 
5 imaging apparatus. A computer 20 was programmed to control 
the operation of the MRI apparatus and process FID signals 
detected therefrom. The gradient field is energized by a 
gradient amplifier 22 and the RF coils for impressing an RF 
magnetic moment at the Larmor frequency is controlled by the 

10 transmitter 24 and the RF coils 26. After the selected nuclei 
have been flipped , the RF coils 2 6 are employed to detect the 
FID signal which is passed through the receiver 28 and then 
through digitizer 3 0 for processing by computer 20. 

Fig. 3 is a schematic representation of a 

15 conventional phased array for use in MRI, as disclosed in 
Roemer U.S. Patent No. 5,086,275. The NMR phased array 10 
includes a plurality of radio frequency (RF) receiver coils 12 
(coils 1 through N c ) defining an imaging volume for the NMR 
imaging of a sample, such as for the NMR medical diagnostic 

20 imaging of a human spine. The separate surface coils 12 are 
identically configured and are arranged in closely spaced 
relationship with overlapping field of view (FOV) , but with 
substantially no interaction between adjacent coils. The 
coils 12 are adapted as part of NMR imaging apparatus such as 

25 shown in Figs. 1 and 2, to simultaneously receive a different 
one of a plurality of NMR response signals, each evoked from 
an associated portion of the sample enclosed in the imaging 
volume. As shown, each coil 12 has its own processing channel 
14 including receiver circuitry 15 and an analog to digital 

30 converter 16. This data processing setup provides for 

constructing a different NMR image for each channel 14 of a 
sample portion from the NMR response signals received by the 
associated coil 12 for that channel 14, and for subsequently 
combining the plurality of different images thus constructed, 

35 on a point by point basis, in the image domain, to produce a 
single, final NMR image of all sample portions from which an 
NMR signal was received by any of the coils 12. As further 
described by Roemer, the optimal combination or weighting of 
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signals from the individual coils 1-N C in the array 10 to 
achieve high signal-to-noise ratio (SNR) is dependent on the 
location of a particular volume element (voxel) . This is 
because the signal of each RF receiving coil C ± is sensitive 
to nuclear spins in proportion to the field B ± created by the 
coil, whereas the noise is "white noise" uniformly distributed 
over the image. Hence, the resultant SNR is a function of 
position. 

As noted above, phase contrast magnetic resonance 
imaging (MRI) uses flow dependent contrast to produce images. 
In the most common method, two data sets are acquired with a 
different gradient first moment in one direction and the pixel 
by pixel phase difference a<£ in the resulting images is used 
to measure the velocity V in the direction of the gradient 
change: a<£ = yaM-lV, where aM x is a changing gradient first 
moment and y is the magnetic ratio. 

In accordance with the invention, a simple and 
efficient method is provided for combining the information 
received from the multicoils that retains the proportionality 
with velocity. 

Consider now phase contrast imaging of a single 
velocity component (two measurements) and let Z xi and Z 2 ± be 
the complex images reconstructed from coil i with sequences 1 
and 2. Similarly, let 

Mi - i\Z xi \ + \Z 2i \)/2 
<J) i = phase (z ljL Z 2l *) 



be the magnitude and phase images from coil i, where 
indicates the complex conjugate and 

phase (Z) - arctan2 (imag(Z) ,real(Z) ) . 
The set {<t> i ) can be thought of as estimates of the true 
velocity induced phase shift a<£. As long as the coils are 
uncoupled the noise in {<p ± } are uncorrelated but have 
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different variances. The variance in 0 ± , c 2 if is inversely 
proportional to K ± 2 : 



Mj/Wf 



The linear combination of {<t> i ) which minimizes the variance 
uses a weight for coil i which is inversely proportional to 
a 2 ± . Thus, one improved algorithm for estimating A<f> is: 



A magnitude weighted image is formed as Ma</> with: 



10 Note that the last equation is the conventional magnitude 
multicoil image. The magnitude weighted phase image is: 



This method requires calculating a separate magnitude and 
phase shift (arc tangent) image for each coil. 

15 An alternate and preferred method combines the data 

from the multiple coils prior to the arc tangent. Each 
complex image, e.g., Z li# contains phase shifts from sources 
other than the desired velocity component. These undesired 
phase shifts are also contained in Z 2 ±, so that the phase in 

20 z n/ z 2i *- s ' ideally, due only to the desired motion. The 

phase of (Z xi Z 2 i*) is exactly the same. For this reason, as 
well as computational efficiency and automatic mapping into a 
2 7T range, the above equation for <t> i is often used to compute 
the phase shift. The quantity (Z lL Z 2 ±) is a unit vector in 

25 the direction <f> ± multiplied by (approximately) K L 2 . Each coil 
provides a similar vector, and, except for noise, all point in 
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the same direction: a<£. Thus, (Z ix Z 2i *) can be viewed as a 
scaled (weighted) estimate of e lA0 . One might therefore ask: 
what is the ideal linear combination of multiple estimates of 
e iA *. To first order, the estimates should have weights 
proportional to M i 2 /W i 2 - Thus we form: 

C = El(z ld z 2i *)/wf] 

M = yf[C\ 

A<p s arctan (imag(C), real(C) ) (lib) 
mag weighted image = M a<£ (11c) 

There have been described methods and apparatus for 
phase contrast MRI with a phased array multicoil in which the 
information from the various coils can be combined while 
retaining the proportionality with velocity. While the 
invention has been described with reference to specific 
embodiments, the description is illustrative of the invention 
and is not to be construed as limiting the invention. Various 
modifications and applications may occur to those skilled in 
the art without departing from the true spirit and scope of 
the invention as defined by the appended claims. 
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WHAT IS CLAIMED IS: 



1. In a multicoil magnetic resonance imaging 
system, a method combining signals detected by multicoils for 
5 phase-contrast imaging comprising the steps of 

a) detecting complex image signals Z u and Z 2i 
with each coil, i, in two excitation sequences with different 
magnetic gradients, 

b) obtaining a combined image C, as the summation 
10 of (Z li Z 2i *)/W i 2 where 

Z 2i * is the complex conjugate of Z 2i 
Wi is a gain factor for coil, i; and 

c) forming an output image dependent on C. 

15 2. The method as defined by claim 1 and further 

including the step: 

d) obtaining phase, <j>, magnitude, M, and magnitude 
weighted image, Ma0, as follows: 

a0 s arctan2 [imag (C),real(c)] 

3. The method as defined by claim 1 and further 
including the step of imaging velocity in a region using said 
combined image . 



20 



25 4. The method as- defined by claim 3 wherein said 

magnetic gradients vary in one direction and said velocity is 
in the direction of gradient change. 



5. The method as defined by claim 1 and further 
30 including the steps of detecting two additional complex image 
signals, Z 3i and Z 4i , and said magnetic gradients vary in 
three dimensions, said complex image signals providing 
components of velocity in said three dimensions. 



35 



6. In apparatus for phase-contrast magnetic 
resonance imaging, the combination comprising: 
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a) means for generating a magnetic field, 

b) means for generating at least one magnetic 
gradient field, 

c) multicoils for detecting complex image signals 
5 Z xi and Z 2i with each coil, i, in two excitation sequences 

with different magnetic gradients; and 

d) means for obtaining a combined image C, as the 
summation of (Z li Z 2 i*) /W± 2 where 

Z 2i * is the complex conjugate of Z 2i 
10 is a gain factor for coil, i; and 

e) means for forming an output image dependent on 

C. 



7. The combination as defined by claim 6 and 
15 further including 

f) means for obtaining phase, 0, magnitude, M, and 
magnitude weighted image, Ma<£, as follows: 

M = /JC[ 

A<p = arctan2 [imag (C),real(C)] 



20 8. The combination as defined by claim 6 and 

further including means for imaging velocity in a region using 
said combined image. 

9. The combination as defined by claim 6 wherein 
25 said magnetic gradients vary in one direction and said 

velocity is in the direction of gradient change. 

10. The combination as defined by claim 6 wherein 
said multicoils detect two additional complex image signals, 
Z 3i and Z 4i , and said magnetic gradients vary in three 
dimensions, said complex image signals providing components of 
velocity in said three dimensions. 



30 

1 



35 



11. In a multicoil magnetic resonance imaging 
system, a method of combining signals detected by multicoils 
for phase-contrast imaging comprising the steps of 
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a) 



10 

detecting complex image signals Z 1± and Z 2i 



with each coil, i, in two excitation sequences with different 
magnetic gradients, 

b) determining velocity induced phase shift image, 
a0, as follows: 



10 



where K± - + \z 2i \)/2 

<t> i = phase (z u Z 2i *) 

= complex conjugate 
Wi = coil gain factor 



12. The method as defined by claim 11 and further 
determining magnitude weighted phase image as follows: 



13. The method as defined by claim 12 and further 
including the step of imaging velocity in a region using said 
magnitude weighted phase image. 

20 14. The method as- defined by claim 11 and further 

including the step of imaging velocity in a region using said 
magnitude multicoil image. 

15. The method as defined by claim 11 wherein said 
25 magnetic gradients vary in one direction and said velocity is 

in the direction of gradient change. 

16. The method as defined by claim 11 and further 
including the steps of detecting two additional complex image 

30 signals, Z 3i and Z 4i , and said magnetic gradients vary in 
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three dimensions, said complex image signals providing 
components of velocity in said three dimensions, 

17. In apparatus for magnetic resonance imaging, 
5 the combination comprising: 

a) means for generating a magnetic field, 

b) means for generating at least one magnetic 
gradient field, 

c) multicoils for detecting complex image signals 
10 Z xi and Z 2 i with each coil, i, in two excitation sequences 

with different magnetic gradients, and 

d) means for determining velocity induced phase 
shift, a0, and M, magnitude multicoil image, as follows: 



15 where Mi - (|z xi | + |z 2i |)/2 

<t> L = phase (Z u Z 2i *) 

= complex conjugate 
Wj_ = coil gain factor 



20 



18. The combination as defined by claim 17 and 
further including means for imaging velocity in a region using 
said magnetic multicoil image, M. 
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